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A Computer-Aided Approach to Nonlinear Control Synthesis 
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*The adaptations accomplished in this project involve inclusion of the describing function capability into version 3.0 of 
the Interactive Controls Analysis (INCA) program. Members of the project are grateful for the technical assistance 
given by Dr. John Downing and Mr. Frank Bauer of NASA-Goddard Space Flight Center who offered help and ideas 
for the framework of the describing function modifications. 



INTRODUCTION AND SUMMARY 

The major objective of this project is to develop a computer-aided approach to nonlinear 
stability analysis and nonlinear control system design. This goal is to be obtained by 
refining the describing function method as a synthesis tool for nonlinear control design. 
This interim report outlines the approach by this study to meet these goals including an 
introduction to the INteractive Controls Analysis (INCA) program which was instrumental 
in meeting these study objectives. A single-input describing function (SIDF) design 
methodology was developed in this study; coupled with the software constructed in this 
study, the results of this project provide a comprehensive tool for design and integration of 
nonlinear control systems. 

Report Outline 

The report consists of several sections in viewgraph format; a table of contents appears 
throughout the report to separate the different sections. In the first section, an introduction 
to the project is found; task information as well as a list of the project accomplishments, 
including an example of INCA's new describing function capability, are given. Next, a 
history of the INCA program is discussed as well as a brief explanation of its capabilities. 
The following section details the implementation of the SIDF analysis capability into 
INCA. The types of nonlinear devices modelled are described and a brief description of the 
theory behind the describing function method is given. The consequent section details the 
nonlinear control synthesis methodology for a flexible spacecraft; steps for developing the 
compensation for a nonlinear control system with flexible modes are presented. This 
methodology is used in the modelling of an actual flight-proven attitude control system, 
INTELSAT 5; this spacecraft employs the pulse-width pulse-frequency (PWPF) modulator 
which is one of the types of nonlinear controllers modelled in INCA. Five different types 
of linear compensation are substituted for the actual compensation; the describing function 
representation as well as the numerical simulation of each sample case is presented. The 



next section details the plan for the second year of the contract; this phase of the contract 
involves much of the same analysis as this past year except the emphasis will be on dual- 
input describing functions. Outlines for possible INCA enhancements are then presented 
followed by a summary and conclusion of this phase of the contract The appendix 
contains supplementary linear analyses of the INTELSAT 5 spacecraft control system 
presented earlier in the report. 

Summary 

The presence of nonlinear devices in a linear system often inhibits linear analysis of that 
system. Nonlinearities discussed in this report come in two forms: energy-storage and 
non-energy storage type. Non-energy storage systems are often representations of 
nonlinear effects which are inherent in the system such as deadzones, limiters, and Schmitt 
Triggers. Describing functions of the energy-storage type include such devices as the 
pulse-width pulse-frequency (PWPF) modulator, the derived-rate modulator, and the 
pulse-width modulator. These types of nonlinear controllers are used to modulate the 
reaction jets of spacecraft attitude control systems and are often referred to as pulse 
modulators. 

The describing function method is often used to predict structural-mode limit cycling which 
may occur as a result of control system/structural interaction. By plotting the describing 
function of a nonlinearity with the corresponding frequency response plot, the designer can 
obtain information about the stability of the system with respect to structural limit cycling. 
For example, the intersection of the two curves indicates the possibility of a nonlinear 
instability; the probability of limit cycle occurring is based on the type of input to the 
system. It is important to note that for energy- storage type nonlinearities, the frequency of 
the intersection point must be the same on both curves. If there is no intersection between 
the curves, approximate phase and gain margins can be defined depending on the relative 




locations of the two arcs. The concept of nonlinear stability margins is important in 
nonlinear control design methodology. 

The describing function capabilities described in this report have been implemented in the INCA 
program. INCA is the result of a software development effort out of NASA Goddard Space 
Flight Center. Using INCA, the designer can execute linear controls analysis in a menu-driven, 
interactive graphics environment; root locus, frequency response, or time response analyses can 
be performed. Furthermore, with the use of a high-resolution graphics terminal, a user can 
query the program for information about the curves, or display multiple plots in a window 
environment Using the modified version of INCA, the controls designer can view describing 
function curves as well as examine the linear/nonlinear interaction in the Nyquist or Nichols 
format for possible limit cycling. Also, INCA's graphics interface facilitates the use of the 
classical gain/phase stabilization concept in nonlinear control synthesis. 

The development of nonlinear control synthesis methodology, especially for the design of on-off 
controller for flexible spacecraft, is perhaps the most important contribution of this project. The 
developed SIDF software was used to demonstrate the practicality of the developed nonlinear 
control synthesis methodology. This procedure was used to compare various linear 
compensators of the INTELSAT 5 yaw-axis RCS controller. 

Step #1: Rigid Body Compensation - The type of compensation for the rigid body 
depends on the performance specifications. This compensation selects the system bandwidth, 
which is related to the rigid body settling time; in general, a lead-lag filter is used where direct- 
rate feedback is not possible. The pole and zero locations (i.e., the lead/lag time constant ratio) 
should be one order of magnitude apart so that the filter does not increase the high frequency gain 


too much. 



Step #2: Pulse Modulator Synthesis - If this item is not fixed as a result of spacecraft 
design, thought must be given to the original configuration of the modulator. Certain factors 
which must be considered in the design such as ease of hardware implementation, stabilization 
requirements of flexible modes, limitations of the reaction control system (i.e. - thruster 
combinations, thrust level). Another consideration is how closely can the nonlinear controller be 
approximated as a linear element. For example, it is difficult to tell by looking at the parameters 
of the PWPF modulator, what its gain and phase characteristics are. Without this knowledge, it 
is difficult for the designer to correlate linear analysis with nonlinear simulation if the modulator 
is assumed to have gain of one and zero-degree phase change. This information can be arrived at 
through describing function analysis of the nonlinearity. 

Step #3: Linear Analysis - Linear Analysis methods such as root locus and frequency 
response will enable the control designer to investigate the stability of the higher-order flexible 
modes. Higher-order compensation may be necessary, in addition to the rigid-body controller, 
to compensate for unstable modes. 

Step #4: Describing Function Analysis - Information about the nonlinearity, as 
approximated in a linear sense, can be obtained through the describing function method. By 
viewing the inverse of the describing function, -l/N(X,co), and the linear plant, G(jco), in the 
frequency domain, information can be obtained about nonlinear stability. For example, if the 
curves intersect, then a high probability for a limit cycle exists; whether or not one occurs 
depends on the type of input into the system as well as other factors. If there is no intersection 
between the curves, approximate phase and gain margins can be defined depending on the 
relative locations of the two arcs. It is important to remember that for energy-storage type 
nonlinearities, the frequency of the intersection point must be the same on both curves. 




Step #5: Nonlinear Simulation - The results in step #4 can be verified by a nonlinear 
simulation, a capability not yet implemented within INCA but which can be accomplished 
through an independent external program. The appearance of a limit cycle can be verified in this 
manner, so can the approximate nonlinear stability margins. Delaying the signal a fixed number 
of cycles throughout the simulation is a simple way of introducing phase lag into the system. 

The phase lag, <|>, follows the relation <J> = 57.3 coT, where co (rad/sec) is the frequency of the 
point of the linear plant close to the describing function and T (sec) is the time delay. By 
varying the time delay, the designer can often induce a limit cycle; the smallest value of <}> which 
causes a limit cycle can be interpreted as the nonlinear phase margin. A similar procedure 
can be followed for the nonlinear gain margin. 

Step #6: Linear Analysis Iteration - If the designer is not satisfied with nonlinear 
stability margins or the possibility of limit cycling, the process can be repeated from step #3 
by altering or adding to the linear compensation. If possible, the controller itself can be 
modified to further decrease the probability of a limit cycle. 

The nonlinear simulation presented in this paper was performed external to INCA. A major 
recommendation of this study is the inclusion of a nonlinear simulation capability within 
INCA. In that way, describing function and nonlinear simulation analyses, both integral 
parts of the nonlinear control synthesis methodology, can be combined in the same 
software. 

A detailed report, which covers these topics in specific detail, is currently under preparation 
and will be available late this year. 
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Summary and conclusions 



The diagram on the facing page illustrates the methodology developed in this study for the design and analysis of 
nonlinear control systems. As can be seen in this figure, the ultimate goal of this project is to develop an integrated 
approach/software for the nonlinear control systems design. 

The user obtains single-input describing function (SIDF) or dual-input describing function (DIDF) plots for a 




Nonlinear Control Analysis/Synthesis/Simulation 



Note: SIDF = Single - Input Describing Function 
DIDF = Dual - Input Describing Function 






In view of the fact that the synthesis of nonlinear control systems is still in the state of development, no attempt 
is made in this project to develop new techniques for the analysis and synthesis of general nonlinear systems. Instead, 
existing techniques such as the describing function method and linear frequency-response methods are employed for the 
development of a computer-aided approach to nonlinear control systems design. 
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* an interim report plus the developed software package 

# final report plus the developed software package 
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This figure represents a generic block diagram for an attitude control system of a spacecraft with structural 
modes. The values given in the viewgraph are for the INTELSAT 5 spacecraft with analog control system (i.e., no 
computational delay , T = 0). However, most spacecraft with digital control system can be represented in this format. 
For example, a similar digital control system has been implemented for the INSAT and ARABSAT spacecraft (see, 
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Pulse-Modulated On-Off Control System 
for a Flexible Spacecraft 



Disturbance 








On the opposite page, a describing function plot of the INTELSAT 5 PWPF modulator is shown. This plot is 
example of the INCA’s describing function capability as a result of modifications made to version 3.0 of the original 
INCA. INCA allows the user to examine several types of nonlinear elements but the software is versatile enough to 
enable interactive analysis of each type of describing function. For example, the opposing figure depicts a describing 
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The opposing plot depicts an example of nonlinear controller interaction with spacecraft's structural modes. The 
figure shows two sets of curves: a PWPF modulator describing function and a four-mode linear plant (i.e., rigid body 
mode plus three flexible modes). It is seen that the two sets of curves intersect at two points. The horizontal lines of the 
describing function represent different input frequencies to the nonlinearity; in order for a limit cycle to exist, the 


V) 

O 

U 
3 S3 
* § 


*■* u 

O 8? 


§• s 

s | 


a JS 
■9 ° 
8 . £ 
c/3 

-s | 

4 — • 

< 03 

4 -* 

O 

c c 

<U 

•a o 

> I- 

•3 a 

8“ H 

J5 1 


C/5 c 

g 5 

b e 

o cx 

5 e 

Q c 

o t; 
c ~ 5 

o y cx 

_ C iu 

C! 3 

O « * 

£ IS 

8 .S 5 

g a .s 


<n w 

■e 


c *"* 
o S g 
^ 'o fe) 

DC- 

■S § -S 

*ts cr TJ 

© e s 

>> <c £ 

O 3 

g © 8 

a & --a 

8^ c o 
2 S3 




PHASE 




Introduction and project overview 
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Summary and conclusions 



INCA is the result of a software development effort out of the Guidance and Control Branch (Code 712) at NASA 
Goddard Space Flight Center. The chief developer of the software is Dr. John Downing while the effort is being directed 
by Frank Bauer. The latest version of INCA, known as version 3.0, is formally distributed through COSMIC meaning 
that it is non-commercial software. INCA allows the designer to perform linear controls analysis in a menu-driven. 
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INCA Data & Command Flow* 



* Adapted from diagram found in INCA 3.0 Primer 










The following four plots represent the graphic interface the INCA user has with the program. These plots are 
taken from the time response and frequency response analyses for linear control systems. The program queries the user 
for the inputs to the systems shown; the values then appear in the block diagram form shown here. This interface allows 
the user to obtain an overall perspective of linear control analysis and design. 
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Impulse . CLOSED LOOP INTER SAMPLE TIME RESPONSE 
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INCA supports describing function analysis of several types of nonlinearities. Many of these arc found in attitude 
control systems; implementation of other common elements is also being planned. The nonlinearities found on the facing 
page are known as "memory-less" systems where the describing functions are not a function of frequency, just amplitude. 
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Schmitt Trigger - The Schmitt trigger is often found in spacecraft attitude control systems using reaction jets, the 
Schmitt trigger logic can be represented as a relay with proportional positive feedback. 




Deadzone Limiter 



The type of nonlinearities presented on the next two pages are also found in spacecraft attitude control systems. 
These nonlinear elements are known as "memory" type since they produce describing functions which are functions of 
frequency as well as amplitude. 


c/3 

a 

<u 

•C 5 

c a 

0 3 

c 2 

1 -i 

.2 ‘2 
wh 3 


O 


3 

a. 

c 


c 

3 

C/3 

C 

o 

o 


3 

a 

C 
• ^ 

<D 


o 

s 

cd 

s 

-C 

o 

Mh 

8- 

Uh 


w 

C 


d« 

S 

13 

c 

o 

i 

cd 

a> 

o 

’€ 

3 

-c 


2 



1 s 


s s 


cd 


w a, 


1 1 


rj 

<S .s 



Nonlinear Control Elements (cont'd) 



Derived-Rate Modulator 










Pulse-Width Modulator - The PWM differs from the previous devices in that it is primarily employed for digital 
controller. Even so, the device's static characteristics can be analyzed similarly to the PWPF modulator and the DRM. 
This device is being considered for the Orbital Maneuvering Vehicle (OMV) Reaction Control System. Like the 
modulators previously discussed, the PWM introduces phase lead into the control loop; the amount of phase lead is 
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Nonlinear Control Elements (cont'd) 



Pulse-Width Modulator 




No attempt is made here to discuss the describing function method in rigor or depth, since this topic is treated 
extensively in many texts (e.g., see Multi-Input Describing Functions and Nonlinear System Design, Gelb 
and Vander Velde, McGraw-Hill, 1968). We will instead emphasize the basic principles and fundamental assumptions, 
and discuss a computer-aided approach to nonlinear control design using the describing function method. 
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Describing Function Method 
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Describing Function Method (Cont'd) 
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The four plots on the facing page show four different states of the PWPF modulator. Shown in the upper left is 
the input signal which corresponds to a sinusoid of to = 5 rad/sec and amplitude, X, equal to .676. This corresponds to 
the leftmost point on the to = 5 curve in the describing function on the previous page. The plot on the lower left 
represents the input to the lag filter, the plot on the upper right depicts the output of that filter which is fed into the Schmitt 
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File: 676. DAT The time is 11:33:48.17 16-JUL-1988 

PWPF Mod. States for X=0.676 

Input Signal After Filter __ 
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Here, the input signal and the first harmonic representation of the output signal is shown. The difference in 
amplitude, gain change, and the change in phase are shown here as well. The gain and phase values define the complex 
value of the describing function N for to = 5 and X = 0.676. The gain and phase values for -1/N are plotted on the 
previous Nichols chart. INCA performs this calculation for the range of user-defined amplitudes and frequencies. 
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The Schmitt trigger and the PWPF modulator are the only nonlinearities represented in INCA that impart phase 
shift to the input sinusoid; as a result, they are easier to view on the INCA screen for they lie off the -180° line. For the 
deadzone, quantizer, and the limiter it is difficult to observe on a black-and-white plot for they lie direcdy on the -180° 
line. 
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The time is 11:01:40.62 4-JUL-1988 

Project : 15 

Schmitt Triggers with +/- uoff 



PHASE 


The describing function depicted here is the yaw-axis PWPF modulator for the INTELSAT 5 (15) spacecraft. The 
plot can be separated into two parts. The roughly vertical line at the right of the plot is known as the minimum-amplitude 
line. This line represents an absolute boundary of the describing function. The series of curves slanting up to the left of 
the graph denote the parts of the describing function determined by altering the input frequency into the nonlinearity. As 
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Summary and conclusions 



The development of nonlinear control synthesis methodology, especially for the design of on-off controller for 
flexible spacecraft, is perhaps the most important contribution of this project. The developed SIDF software was used to 
demonstrate the practicality of the developed nonlinear control synthesis methodology. This procedure was used to 
compare various linear compensators of the INTELSAT 5 yaw-axis RCS controller. 
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control designer to investigate the stability of the higher-order flexible modes. Higher-order compensation may be 
necessary, in addition to the rigid-body controller, to compensate for unstable modes. 
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Step #4: Describing Function Analysis - Information about the nonlinearity, as approximated in a linear sense, can 
be obtained through the describing function method. By viewing the inverse of the describing function, -l/N(X,(o), and 
the linear plant, G(jto), in the frequency domain, information can be obtained about nonlinear stability. For example, if 
the curves intersect, then a high probability for a limit cycle exists; whether or not one occurs depends on the type of input 
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Step #6: Linear Analysis Iteration 



The INTELSAT 5 yaw-axis RCS controller was used to validate the SIDF analysis methodology developed in this 
project. The information for the plant, linear compensation, and nonlinear controller is listed on the opposing page. The 
linear compensator parameters listed in the block diagram are for case #1. 



Pulse-Modulated On-Off Control System 
for a Flexible Spacecraft 



Disturbance 



Five different linear compensators are summarized on the facing page. In cases #1 and 2, simple lead-lag 
compensation is employed to stabilize both rigid body and flexible modes. Cases #3,4,5 utilized different structural 
filtering to avoid the structural mode limit cycle predicted by the describing function method. 



INTELSAT-5 RCS Design Examples 
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- Case #1 with z=3.4, t, z =0.002, p=3.4, £, p = 1 

Case #5: Lead-Lag Compensator with Notch Filtering 

- Case #1 with z=3.6, C z =0.002, p =3.6, £ P =1 







The left plot displayed on the opposing page shows the small nonlinear gain and phase margins of case #1 with 
respect to strucural mode limit cycling. The plot on the top right is close-up view of the area in the vicinity of the 
describing function. The plot on the bottom right side of the page is also a close-up view of case #1 with 100 msec time 
delay; structural mode limit cycling at 3.5 rad/sec is predicted from describing function analysis. 
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These plots represent the nonlinear simulation of case #1 without delay; note the absence of limit cycling in the 
plot on the lower right. 
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These plots represent the nonlinear simulation of case #1 with delay; note the appearance of limit cycling in the 
plot on the lower right. 
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Introduction and project overview 
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Summary and conclusions 



Although the SIDF method is a useful tool for certain nonlinear systems, the so-called DIDF (Dual-Input 
Describing Function) method is needed in certain cases where two sinusoidal signals appear at the input of a nonlinear 
element. In these cases, the high frequency signal is taken to be the system limit cycle, while the slowly varying signal i: 
coming from input commands or disturbances to the system. The concept of the DIDF is rather straightforward, but the 
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The plots on the opposing page represent a variation of the cases presented earlier in the report. The example 
presented here depicts the INTELSAT 5 spacecraft with the lead-lag filter used in case #1; in this example, a constant 
torque of 0.5 N-m. is introduced into the system. Note the steady-state value of yaw angle with the "beating” effect. 
More importantly, notice the lack of limit cycling in the torque output plot. The thruster is counteracting the 0.5 N-m 
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This system is nearly identical with the one of the previous page except that a torque bias of 0.2 N-m is introduced 
as indicated by the lower duty cycle firing of 10%. The torque command has a slightly negative average bias; as a result, 
some positive thrusting occurs. 



Yaw An g 1 e (deg) 
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This system again represents case #1 except that torque and attitude commands are input to the system. As a 
result, the steady-state torque output looks much as the previous system did. 
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Future INCA Enhancements 
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